Electronic structure and transport in CsBi4Te6. 
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CsBi4Te6, is calculated and analyzed using the semi-classic transport equations. It is 
shown that to obtain a quantitative agreement with measured transport properties 



\j ' a band gap of 0.08 eV must be enforced. A gap in reasonable agreement with exper- 

iment was obtained using the generalized gradient functional of Engel and Vosko. 



We found that the experimental p-type sample has a carrier concentration close to 
optimal. Furthermore the conduction bands have a form equally well suited for 
thermoelectric properties and we predict that an optimally doped n-type compound 



T3 . 
rH , could have thermoelectric properties exceeding those of the p-type. 
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O ■ I- INTRODUCTION 

O ' The search for new thermoelectric materials is a quest to maximize the dimensionless 

^ . 

^ ' figure of merit zT = [aT / k)S'^^ where S is the Seebeck coefficient and a and k are the 

' electronic and thermal conductivities respectively. zT quantifies the performance of a ther- 

moelectric and one must therefore maximize the power-factor S^a and minimize k. As S, 
a and k are coupled and all depend strongly on the detailed electronic structure, carrier 
concentration and crystal structure, finding new compounds with large values of zT is a dif- 
k> ■ ficult task. Still, several completely new types of materials, with complex crystal structures 

j^ . and 2;T's exceeding the presently used alloys, have recently been found.^^^ 

CsBi4Te6 is remarkable because of its thermoelectric properties at low temperatures 
{^zT ^ 0.8 at 225 K).^ It has a complex layered crystal structure and some intriguing direct 
Bi— Bi bonds.— It would therefore be interesting if a direct link between the bonding and 
band structure of CsBi4Te6 and the thermoelectric properties could be estabhshed, as has 
been done earlier for other thermoelectric compounds . ^i^i^i^i^ One study of the band struc- 



FIG. 1: (color online) Structure of CsBi4Te6. The structure consists of long parallel rows of 
BiTeg octahedra (green/light grey) stitched together by TesBi— BiTes bonds (blue/dark square 
based pyramids and bonds). The BiTe-layers are separated by layers of Cs atoms which act as 
electron donors to the framework. The unit cell shown is the C2/m monoclinic cell a = 51.9205 A, 
h = 4.4025 A, c = 14.5118 A and /3 = 101.4800°. 

ture of CsBi4Te6 has been publishedi^ and very recently this study was used to interpret 
the resuhs of a systematic experimental study of CsBi4Te6 doping.— Though the earlier 
band structure study did report the effective masses^", we wish to improve the quantitative 
link between the band structure and the thermoelectric quantities. Furthermore, despite 
using the same LAPW based method, we obtain a somewhat different band structure than 
was found earlier.^^ We suspect that this is due to relativistic effects being poorly treated 
in the earlier study.— These are extremely important due to the large spin-orbit splitting 
of the Bi-6p valence states. The challenge of obtaining a correct band structure for the 
bismuth-telluride compounds is probably best illustrated by the attempts to calculate the 
band structure of the well-known thermoelectric compound Bi2Te3. De-Hass-van Alphen 
experiments^^ and angle-resolved photo-emission studies^^ have shown that both the lowest 
conduction band (LCB) and the highest valence band (HVB) have six-fold degenerate band 
edges. Only recently has this been reproduced by band structure calculation o^i^^i^^ and only 
very recently has also the gap been correctly predicted.— 

The paper is organized as follows: First we briefly review the crystal structure and 
measured thermoelectric properties. We then discuss the size of the band gap and compare 
the theoretical prediction with the experimental measurements. The dependency of the 
transport properties on the carrier concentration is then reported and finally we discuss the 
band structure and bonding. 

II. STRUCTURE AND MEASURED TRANSPORT 

The structure of CsBi4Te6 is shown in Fig. ^ The atoms are arranged in Cs layers and 
Bi4Te6 slabs with Bi— Bi bonds stitching the slabs together along the a-axis. The Bi— Bi bond 
is unusual as it involves a reduction of Bi'^^ to Bi^+ seldom seen in Bi-chalcogenide systems.-^ 
The inner located Bi atoms are coordinated to six Te atoms in a distorted octahedral whereas 



Bi atoms involved in Bi— Bi bonding are coordinated to five Te atoms. 

CsBi4Te6 has a side centered monoclinic unit cell. In the C-centered setting the a-axis 
is considerably longer than the b- and c-axes, Fig. ^ The slabs stretch out infinitely along 
the 6-axis making CsBi4Te6 a ID-like structure along the short 6-axis. Von Neumanns 
principle relates the point group symmetry of the crystal structure (monoclinic) to the 
tensor properties of the crystal. The monoclinic symmetry with f3 ^ 90° means that the 
conductivity tensor will have the symmetry 
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and similarly for the Seebeck coefficient. 

Experimentally a good conductivity was only found parallel to the short 6-axis,iii as 
would be expected from the crystal structure. Consequently good thermoelectric properties 
were only found when measured parallel to the 6-axis.^'^^ The measured Seebeck coefficient 
rises monotonically up to approximately 275 K (S* ~ 175yuV/K) and starts to decrease at 
higher temperatures. This is in good agreement with the very narrow band gaps, ranging 
from 0.04 to 0.1 eV, found experimentallyi»ii*i^ and means that CsBi4Te6 only has favorable 
thermoelectric properties at low temperatures. It is also important for the thermoelectric 
properties of CsBi4Te6 that a high electrical conductivity (~ 1450 S/cm at 250 K) and 
an extremely low lattice thermal conductivity of k/ = 0.6 W/mK (calculated by using 
the Wiedemann-Franz law to subtract He from the measured n) parallel to the 6-axis were 

III. COMPUTATIONAL METHOD 

A. Electronic structure 

The calculations were performed using the L/APW+lo method^- as implemented in the 
WIEN2k code.^*^ A plane wave cutoff defined by min(i?a)max(fc„) = 6 and sphere sizes 
of 2.7 were used. The exchange-correlation potential was calculated with the Perdew- 
Burke-Ernzerhof (PBE)^^ and the Engel-Vosko (EV)^ generalized gradient approximations 
(GGAs). 36 /c-points on a shifted mesh in the IBZ were used for the SCF calculations 



As will be discussed later we find a band gap of zero using the PBE-GGA. This disagrees 
with the earlier band structure study of CsBi4Te6, which found a gap of 0.04 eV.'-^ A gap 
of 0.04 eV is in reasonable agreement with experiment, but we believe this was caused by a 
cancellation of errors. We can see two possible sources of error in the previous paper: either 
a poor treatment of spin-orbit coupling or an inadequate fc-point sampling. 

In WIEN2k spin orbit coupling is included in a second variational step.^^ The size of the 
second variational basis set is controlled by an energy cut-off which limits the number of 
eigen vectors used. Furthermore pi/2 local orbitals can be included in the second variational 
step.— Because of a strong spin-orbit splitting of the Bi 6p-states at the Fermi-level, a correct 
description of relativistic effects is extremely important and the second variational basis-set 
must be well converged. We first performed a calculation using the default energy cut-off 
(1.5 Ry) of the WIEN code and did find that this induces a small gap of 0.012 eV. To test 
convergence we then performed two calculations: one setting a high energy cut-off (6.0 Ry) 
and another with a smaller cut-off (2.0 Ry) but including pi/2 local orbitals^^ for the Bi 
atoms at energies close to the Fermi level. Both calculations found a band gap of zero 
and a very good agreement between the two calculations was found. No details about the 
second variational basis set were reported in the earlier paper,»12- but it should be pointed 
out that a published study on Bi2Te3,— had to be corrected due poor treatment of the pi/2 
contribution.— 

The earlier papei»^ found the LCB minimum at the (0.881,0.881,0.175) position in the 
primitive reciprocal basis used in that study. This is the (0.238,0,-0.175) position in the 
C-side reciprocal basis used here. The earlier minimum position is not in complete agree- 
ment with the minimum position found here: (0.259,0,-0,166). The minimum of the LCB 
calculated band gap is not at a high symmetry position which means that an adequately 
dense fc-mesh is very important for the exact location. The earlier paper does not report 
the procedure used to find the minimum, but only reports using 13 k points in the IBZ for 
the self-consistent calculation, which is certainly to coarse a mesh for an exact location of 
the LCB minimum. We determined the minimum on the 1330 k point mesh used in the 
transport calculations discussed below. 
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FIG. 2: Integrand factors df/de and (e — ij)df /de in Eqs. Q-® at T = 250 K. Arbitrary y units. 

B. Transport 

The band structure was analyzed using semi-classic Boltzmann theorj*24 and the rigid 
band approach. The rigid band approach to conductivity is based on the transport distri- 
bution 

^"/^(^) = ^E^"/^(^'k)^^^^ (2) 
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The k-dependent transport tensor is given as 



crc,/3(i, k) = eVi,kt;„(z,k)t;^(i,k) 



(3) 



where r is the relaxation time and ^^(zjk) is a component of the group velocities. The 
transport coefficients can be calculated as a function of temperature and chemical potential 
by integrating the transport distribution 
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The bands, and hence (y{e), are left fixed (thus "the rigid band approach") and therefore 
only one band structure calculation needs to be performed per compound. The number 
of carriers is changed by varying the chemical potential in Eqs. (jHEj). Fig. |21 shows the 
integrand factors in Eqs. (jHEI) at T = 250 K. It is seen that the distribution is quite broad 
and the transport coefficients are thus a sum over several Fermi surfaces. It is therefore very 
important that the band is correctly calculated. 
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FIG. 3: The bb component of the Seebeck coefficient tensor at various fixed band gaps, (a) Sbb vs 
carrier concentration at a fixed temperature of T = 250 K. The dotted horizontal hne marks the 
Sbb = 175 fj-Y/K found experimentally at 250 K. (b) Sbb vs temperature. The carrier concentration 
for both curves has been fixed so that S = 175 fiY/K at 250 K. 

In Eqs. (jHl) the relaxation time is unknown. Our approach in the present work is too treat 
it as a constant. The Seebeck coefficient, S = a'^u, is then independent of r and can thus 
be calculated on an absolute scale. The conductivity can only be calculated with respect to 
the relaxation time. As ni and r are not readily available from band structure calculations 
they must somehow be included as parameters. 

For the transport calculation, eigen energies at 1330 A;-points on a non-shifted mesh in 
the IBZ were calculated. For calculation of the necessary derivatives, Eq. (jH)), the program 
BoltzTraP was used.'^ BoltzTraP relies on a well tested smoothed Fourier interpolation to 
obtain an analytical expression of the bands. ^^ The original fc-mesh was interpolated onto a 
mesh five times denser then the original. 

For computational reasons our calculations were carried out in the 5-centered setting of 
the cell with a = 51.9205 A, 6 = 51.0530 A, c = 4.4025 A and 7 = 163.8256°, where the 
c-axis in the S-centered setting is parallel to the 6-axis in the C-centered cell. However, we 
have converted all transport properties and band structures back into the basis of the C- 
centered cell and will discuss the transport properties parallel to the 6-axis of the C-centered 
cell. 



IV. RESULTS AND DISCUSSION 



A. Band gap 



The calculated band structure of CsBi4Te6, which will be discussed in detail later, is pre- 
dicted to have an indirect band gap of zero. This disagrees with the band gap of 0.04—0.1 eV, 
found experimentally.— »ii*i^ As a result of the underestimated band gap the calculated See- 
beck coefficient of the standard PBE-GGA calculation is too low at elevated temperatures. 
This is illustrated in Fig. IHK, where the calculated S (full line) does not agree with experi- 
ment, irrespective of carrier concentration. 

The underestimation of the band gap is a well known problem of Kohn-Sham theory . ^^i^'i'i^^ 
Experience suggests that the shape of the bands are correct and that most of the problem 
can be amended by a rigid shift of the conduction band. We therefore introduced a gap 
between the minimum energy of the LCB and the highest energy of the HVB by hand. The 
magnitude off the gap was then estimated by comparing the calculated Seebeck coefficient 
with the experimental values.^ The first observation needed to fix the gap was that a gap 
larger than 0.06 eV is needed in order to obtain the experimental value oi S = 175 /xV/K 
at 250 K, Fig. Et^. Secondly, experiments have shown that the Seebeck coefficient starts 
to decrease at temperatures above approximately 250 K.^ Fig. Eb shows that with a gap 
of 0.09 eV the Seebeck coefficient continues to increase up to approximately 300 K which 
gives an upper-bound of the gap. With a gap of 0.08 eV the experimental behavior is well 
reproduced and in good agreement with the band gap of 0.04-0.08 eV estimated from the 
experimental transport coefficients.— 

B. The Engel-Vosko gap 

The problem of underestimating the band gap in DFT is often discussed and improve- 
ments have been suggested using the GW approximation2L22iS2i^ or by including exact 
exchange.— >^i*^2i^ In some previous papers^'^^ it has been suggested that one route to a bet- 
ter band gap could be the EV-GGA.^" It is known to give very poor total energy differences, 
but for a number of narrow band gap semiconductors the calculated band gaps (i.e 1.23 eV 
for Si, 0.39 eV for Ge and 1.04 eV for GaAs) compare much better with experiment than 
the PBE-GGA. As a recent example we predicted a band gap of 0.4 eV for the type- VIII 
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FIG. 4: (a) Power factor S'^a vs doping at 250 K. (b) Figure of merit (zT) vs doping at 250 K. 
The parameters r = 13.7 x 10~^^ s and ki = 0.6 W/mK were used (see text). Ke was calculated 
using the Wiedemann-Franz law. The vertical line indicates the carrier concentration that agrees 
with experiment. 

EugGaigGeso clathrate^ which was subsequently confirmed experimentally.— 

We have applied the EV-GGA to CsBi4Te6, retaining the other computational parameters 
as reported in section ITTTl resulting in Eg = 0.06 eV, in good agreement with both exper- 
iment and the above analysis. To illustrate the usefulness of the EV-GGA we have also 
performed calculations on Bi2Te3. Using the PBE-GGA we obtain a band gap of 0.11 eV in 
agreement with an earlier calculation.— Using the EV-GGA, but otherwise exactly the same 
computational parameters as earlier—, we obtained a gap of 156 meV in excellent agreement 
with Eg = 154 meV calculated using the sX-LDA^^ functional and in good agreement with 
the zero-temperature extrapolated experimental value of 162 meV.— 

The consistency of the better band gaps with the EV-GGA is thought provoking and 
could hint to an underlying reason. A possible explanation could lie in the construction 
of the EV-GGA which put used the virial relation is used to construct a functional which 
reproduces atomic exchange potentials better than the usual functionals.— It is known that 
the one-electron potential in the atomic limit jumps disco ntinuously at integer electron 
values'^^, and that this discontinuity can give a large contribution to the band gap.— It is 
possible that the EV-GGA, where , also gives a better reproduction of this discontinuity. 
However, the better band gaps of the EV-GGA are a phenomenological observation that 
should warrant further investigation. 



C. Optimal carrier concentration 

Effective masses (m*) can be calculated from the band structure as a function of car- 
rier concentration and temperature.-^ The effective masses are related to the conductivity 
through the relaxation time approximation: a = e^{n/m*)T. From the band structure we 
obtain e^n/ml^ = cXbb/T = 1-06 x 10^^ {Q m s)~^ at the estimated experimental carrier 
concentration. Using the experimental conductivity of abb = 1450 S/cm at 250 K we ob- 
tain a relaxation time of r = 13.7 x 10~^^ s. Thus o"b(,/r is not unusually high but r is 
substantially larger than what has been found for Bi2Te3^ and comparable to r in simple 
metals. Keeping r constant for all doping levels, we calculate the power factor as a function 
of doping, Fig. |3^. It can be seen that CsBi4Te6 has a high power factor comparable to that 
of Bi2Te3 and that the estimated experimental carrier concentration is close to optimal for 
a p-doped sample. Fig. ^p shows the calculated zT using the experimental lattice thermal 
conductivity, ki = 0.6 W/mK. 

Fig. E^ also reveals that the optimal power factor and zT for n-doped CsBi4Te6 is even 
higher than for optimally p-doped. The high estimated power factor and zT for n-doped 
CsBi4Te6 is interesting considering that n-doped CsBi4Te6 samples have been successfully 
synthesizediiii and a thermoelectric device needs both p- and n-type legs. Experimentally 
the ra-type compounds are found to have a smaller absolute Seebeck coefficient and a higher 
conductivity. This can be expected if the carrier concentration is higher than in p-type 
compounds. However, this explanation disagrees with the observation that the Seebeck co- 
efficient starts to decrease at a lower temperature for the n-type than for the p-type. This 
can only be explained if the n-type materials have a larger concentration of impurities, giv- 
ing rise to states in the band gap. Furthermore it should be pointed out that electrons and 
holes couple differently to the lattice vibrations. At 250 K the scattering mechanisms are 
dominated by electron-phonon coupling and one should use different relaxation times for 
electrons and holes. Unfortunately the poor agreement between the calculated and exper- 
imental temperature dependence of S makes it difficult to estimate a carrier concentration 
for the n-doped samples. The results for n-doped materials should therefore be taken with 
some caution but the results do suggest that further attempts at optimizing the doping for 
n-type materials could be worthwhile. 
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FIG. 5: Band structure of CsBi4Te6. The points are labeled with respect to 7r(a*,6*,c*) where 
a* , }f and c* are the reciprocal lattice vectors of the C-centered unit cell (see Fig. EJ . The band 
gap fixed at 0.08 eV. 

FIG. 6: Constant energy surfaces: (a) 0.05 eV below the optimal chemical potential of the HVB 
(see Fig. EJ. (b) 0.05 eV above the optimal chemical potential of the LCB (see Fig. inj. The lattice 
vectors of the reciprocal C-centered unit cell are shown. 



D. Bonding and band structure 

Fig. El shows the band structure of CsBi4Te6 with Eg = 0.08 eV. It is seen that the 
maximum of the HVB is found at the F-point while the minima of the LCB he in the 
a*c* plane at the (0.259,0,-0.166) and (-0.259,0,0.166) points. With the indirect gap set to 
0.08 eV, a direct band gap of 0.1 eV is found, in good agreement with the optical absorption 
measurements which also found a gap around 0.1 eV.— 

It can be shown by a simple argument, that the power factor grows with increasing slope 
and decreasing height of the transport distributioii^ and good thermoelectric properties can 
be related to a large anisotropy in the transport distribution around the chemical potential. 
This is caused by the anti-symmetric shape of {6~fi)df /de, Fig.|2l that enters the expression 
for the Seebeck coefficient. The optimal carrier concentration will therefore correspond to 
a chemical potential close to the band edges. The extrema of the {e — ^)df /de function 
at 250 K are located approximately at je — /i| = 0.05 eV, Fig. |21 From the shape of 
{e — ^)df /de two favorable features for thermoelectric performance can be inferred: (i) 
The constant energy surface approximately 0.05 eV from the band edge should have a large 
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(a) (b) 
FIG. 7: (color online) Constant electron density surfaces calculated from (a) the highest valence 
band state at the T-point and (b) the lowest conduction band at the band edge (the (0.260,0,-0.166) 
point, see text). The surface covers a volume where the electron density is larger than 0.005 e/A. 
The small spheres in the middle are the Cs, the dark grey /blue sphere represent Bi and the light 
grey /green spheres represent Te. 

area which decreases rapidly when approaching the band edge and (ii) the carriers at these 
energies should have reasonably high group velocities to maximize the transport distribution, 
Eqs. (j21ini)- Normally, the area of the Fermi surface is inversely proportional to the group 
velocity, so the two conditions above can only be achieved if there is a large anisotropy in 
the band structure or several carrier pockets contributing to the conductivity. 

From Fig. it is clear that the large surface area is due to the small dispersion of the 
LCB in the a*c*-plane, especially along the a* axis, and the HVB along the a* direction. 
This is also illustrated in Fig. IHl which shows the constant energy surfaces at 0.05 eV below 
and above the HVB and LCB edges. These low dispersion directions co-exist with large 
dispersions along the h* axis leading to a good conductivity parallel to the 6-axis in direct 
space. 

The chemical question is then obviously why the dispersion along the a*-axis is small. To 
analyze this the constant electron density surface was calculated from the HVB state at the 
F-point, Fig. [7| As was also pointed out earlier— the density at this state is mainly situated 
in sheets along the a-axes and it turns out that the density only changes very little for the 
HVB state at a*, resulting in the small energy dispersion. As opposed to what was reported 
earlier, also the minimum of the LCB is associated with the same atoms, but now with a 
clear anti-bonding character. 

V. CONCLUSION 

We have calculated the electronic structure of the promising thermoelectric compound 
CsBi4Te6 and shown that a band gap of 0.08 eV must be enforced to obtain a quantitative 
agreement between the calculated and measured transport properties. It was also shown 
that a gap in reasonable agreement with experiment can be obtained using the EV-GGA 
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functional. We have found that the experimental p-tjpe sample has a carrier concentration 
close to optimal. The conduction bands have a form equally well suited for thermoelectric 
properties and we predict that the n-type compound might have thermoelectric properties 
exceeding those of the p-type. 
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